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Abstract 

We investigate scaling and correlations of the energy and momentum in an 
evolving network of cosmic strings in Minkowski space. These quantities are of 
great interest, as they must be understood before accurate predictions for the power 
spectra of the perturbations in the matter and radiation in the early Universe can 
be made. We argue that Minkowski space provides a reasonable approximation to 
a Friedmann background for string dynamics and we use our results to construct a 
simple model of the network, in which it is considered to consist of randomly placed 
segments moving with random velocities. This model works well in accounting for 
features of the two-time correlation functions, and even better for the power spectra. 
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1 Introduction 



Phase-ordering dynamics after a cosmological phase transition may explain the origin of 
structure in our universe [0J . Our interest is in the role of one class of theories - cosmic 
strings 0. 

Recent attention has focused on cosmic strings formed during the breaking of the 
grand unified symmetry, which naturally generates perturbations in the matter and ra- 
diation fields of the right order for structure formation. However, to make diagnostic 
predictions between cosmic strings, other defects and inflation-based theories we need 
to use a detailed theory of the cosmological evolution of perturbations. There exists a 
well-founded theory for the standard inflationary scenario of gaussian fluctuations in the 
gravitational potential. In extending the theory to cover defect-based theories there are 
three principal issues to be faced. Firstly, although the initial conditions are random, 
they are not gaussian which makes averaging over the ensembles of defects a numerical 
problem. Secondly, the creation of the defects must conserve energy-momentum: hence 
there are compensating fluctuations in the fluid components on super-horizon scales. 
Thirdly, in the subsequent evolution, the defect stress-energy 0^ must be included in 
the Einstein equations. Defects form an unusual component of the total stress-energy, 
in that the metric perturbations only affect the evolution of the defect network to sec- 
ond order, so in the linear approximation we can evolve the network in the unperturbed 
background metric. Its stress-energy then acts as a source term for the familiar set of 
differential equations governing the evolution of the perturbations. As demonstrated by 
Veeraraghavan and Stebbins ||, the solutions to the Einstein equations describing the 
subsequent evolution of perturbation variables with a source term can be expressed as a 
convolution of a suitable Green's function (dependent on the usual cosmological param- 
eters) with the source term integrated over time. For a set of perturbation variables A a 
and a set of source terms Sb related to Qa U 



where r\ and rf are conformal times and k is the wavenumber of the mode. Calculating 
the power in these variables will then involve integrating over the two-time correlator 



When using the Veeraraghavan- Stebbins formalism, the Green's functions can be calcu- 
lated from standard cosmological perturbation theory: analytically in the simplest cases, 
or numerically for a more realistic universe. However, in the absence of a workable ana- 





(|A a (M)| 2 > = p drj' drj" G ab (r], r]')G ac (rj, rj") (S b (k, r}')S c (k, r]")) 



(2) 
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lytic treatment of string evolution (although see |3J] and ||) we must rely for the moment 
on a numerical study of the two-time correlators. 

In a simplified model the Universe consists of two fluids: a pressureless CDM com- 
ponent and a tightly-coupled baryon-photon fluid. In the convenient synchronous gauge, 
the sources are Goo an d ikjQ j (which we shall hereafter call U), when the equations 
are rewritten using the energy-momentum pseudotensor introduced by Veeraraghavan 
and Stebbins |£§. In the Newtonian gauge they are also the sources of curvature fluctu- 
ations ||. It is therefore the two-time correlation functions of these variables that we 
study. Such a study is particularly timely in view of recent work on the power spectrum 
of Cosmic Microwave Background around the degree scale from strings || and from 



global textures [g, [L0[]. As emphasised by Albrecht et al |nj, the distinctive appearance 
of "Doppler" peaks and troughs seen in inflationary calculations and texture models de- 
pend sensitively on the temporal coherence of the sources. If one assumes little coherence 
the peaks are washed out; on the other hand an assumption of total coherence preserves 
them. Magueijo et al in [0] assumed that strings were effectively incoherent and obtained 
a rather featureless CMB power spectrum at large multipole I. In || this assumption 
was justified by a numerical study of the two-time energy-density correlator. 

In contrast both Durrer et al || and Crittenden and Turok |l(| assumed maximum 
coherence for their texture models and found that the peaks were preserved, albeit in 
positions typical of isocurvature models 0. It is therefore clear that understanding 
the temporal coherence of string sources is of great importance when calculating their 
microwave background signals. 

In this paper we present the results of some numerical "experiments" evolving strings 
in Minkowski space. Using Minkowski space rather than a Friedmann model may seem 
rather unrealistic. However, we expect a network of strings evolving in Minkowski space 
to have all the essential features of one in a Friedmann background. A tangled initial 
state consisting of a few probably "infinite" strings plus a scale free distribution of loops 
straightens out under tension and continually self intersects resulting in the transfer 
of energy into very small fast moving loops. The infinite string approaches "scaling" in 
which it is characterised as a set of more or less Brownian random walks. The step length 
of the walks and their average separation are both approximately equal to an overall 
network scale £, which increases linearly with time (as naive dimensional analysis would 
suggest). The length density of the infinite string decreases as £ -2 , again as dimensional 
analysis suggests. In fact, £ is conventionally defined so that the length density is precisely 
l/£ 2 The effect of an expanding background is to damp coherent motions on super- 
horizon scales. However the network scale is much less that the horizon size so one can 
argue that the expansion does not significantly affect the network dynamics. The great 
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advantage of Minkowski space is that the network evolution is very easy to simulate 
numerically: the code is generally many times shorter than an equivalent code for a 
Friedmann background [[13], [14|], and makes fewer demands on both RAM and CPU time. 
We are therefore able to get much better statistical significance from the data. 

We present results from an extensive programme of numerical simulations for the two- 
time correlators of 0oo and U, and use a simple model to explain most of the features we 
observe. It turns out that, to a good approximation, the network can be thought of as 
consisting of randomly placed segments of string with random velocities. The segments 
are of length £ and number density £~ 3 , thus reproducing the correct behaviour for the 
length density. 

In this model it is easy to show that the coherence time-scale in a Fourier mode of 
wavenumber k is determined by the time segments take to travel a distance The 
model in fact predicts that the correlations between the energy density at times r\ and 
rf decrease as exp(— ^k 2 v 2 (rj — r/') 2 ), where v is the r.m.s string velocity. Hence we can 
talk of a coherence time-scale r/ c = \/3/vk, which is the time over which the correlation 
function falls to of its equal time value. Given that we find v 2 = 0.36, the model 
predicts r] c ~ 3/k. Our results actually indicate that at high k, r] c decreases faster than 
k~ l , but we have some evidence that this behaviour is a lattice artifact. 



2 Flat Space Strings 



We use a development of a code written by one of us some time ago |T5|, [L6[]. Initial string 
realisations are generated using the Vachaspati-Vilenkin algorithm |]17|| . This mimics 
the Kibble mechanism for the spontaneous breaking of a scalar U(l) symmetry with a 
'Mexican hat" potential, as the Universe cools through a thermal phase transition. 

This initial configuration is set up for evolution on a cubic lattice with fundamental 
lattice side 5. We are free to alter the initial correlation length £0 i n terms of S and we 
may add structure to the network by placing a percentage P c of cusps randomly along 
the string. Cusps are string links which are confined to one lattice site and which move 
at the speed of light. Adding cusps avoids peculiarities arising from having long straight 
segments of string in the initial conditions. 

In Minkowski space the string equations of motion are 

X" = X (3) 

where ' = d/drj and ' = djdo\ rj is time and a is a space-like parameter along the string. 
If we regard Minkowski space as a Friedmann-Roberston- Walker Universe in the limit 
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that the expansion rate goes to zero, then rj can be identified as FRW conformal time. 
X = X(<7, if) is a position three- vector which satisfies 

X' ■ X = (4) 

X /2 + X 2 = 1 (5) 

This is a convenient set of constraints, as ((5|) ensures that the energy of a segment of 
string is proportional to its length in a. 

We solve this using the Smith- Vilenkin algorithm |TJ| which is based on the exact 
finite difference solution to (|3|), 

X(a,r) + 5) = X(a + 6,T])+X(a-5,r))-X(a,r)-5) (6) 

If the string points are initially defined on the sites of a cubic lattice (N5) 3 , then (^) 
ensures that they remain on the lattice at time steps of 5. One can see from (Q) that 
stationary elementary segments have length (and hence energy) of 25 and point in one 
of six directions. 

Because the string points lie on the sites of the lattice, identifying crossing events 
is easy. When two strings cross, they intercommute with a probability which is set to 
1. Loops with energy greater than or equal to a threshold value of E c are allowed to 
leave the network, while reconnections are forbidden for loops with energy equal to E c . 
Forbidding reconnections allows energy to leave the network fairly efficiently; otherwise 
it takes much longer for the effect of the initial conditions to wear off. The natural and 
usual value for E c is the minimum segment length, 25. In fact, most of the energy in the 
string network goes into the smallest possible loops. In some sense, the cusps model the 
gravitational radiation of a realistic string network: They travel at the speed of light and 
do not subsequently interact with the network. 

Each realisation is evolved on a (645) 3 or a (1285) 3 lattice with approximately 7500 or 
60000 points describing the string network. For calculating quantities like power spectra 
and two-time correlation functions we typically average over 50 realisations. 



3 Results 

The energy-momentum tensor for a cosmic string at a point x is 

6^(x) = J da(X^X u - X'^)5 (3) (x - X(a,r))) (7) 

It is simple to calculate this on a cubic lattice for our string realisations, and use a Fast 
Fourier Transform to get a Fourier mode representation. Here we want to measure the 
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Figure 1: The energy-density length scale £ over time rj for a 128 3 lattice, averaged over 
50 realisations 



power spectra of Qoo{k) and U(k). These are calculated by averaging the amplitudes 
over all Fourier modes with a wavenumber between k and k + 2%/ 8. 

Previous work [13, R, [Tj| has tried to identify those features of the network that are 
scaling, when a relevant length scale grows with the horizon. For example, we can define 
the familiar energy density scale £ by £ 2 = /i/pmf where pi n f is the density of string with 
energy greater than £. (This apparently circular definition works because we know the 
initial step length £ and to calculate p in f we use £ from the previous time step). Scaling 
is reached when x = d^/dri is constant. In the cases considered here x = 0.15 (E c = 28) 
and x = 0.12 (E c = 48). The usual definition for scaling, that x = £/?7 is constant, can 
not be used here as throughout our simulations £o is too large to be ignored. For this 
reason, we express all scaling functions in terms of £, rather than rj. The behaviour of £ 
is shown in Figure [l]. 

The infinite string energy-momentum power spectra also exhibit scaling behaviour, 
which we express in terms of the scaling functions P p and P u . We also consider the 
equal time cross correlation function (U(k)Ql (k)) which is a real function within the 
ensemble errors. 

<|eoo(M)| 2 > = ^^ (8) 
(u(k)e* (k)) = VXPU ^ (9) 

(\U(k^) = ^^l (10) 



6 




Figure 2: Scaling function for the energy density as defined by equation The solid 
line is the measured scaling function with 1 — a error bars from ensemble averaging. The 
dashed line is the fitted function equation fll2|). The dash-dot line is the predicted form 
from our model, the last line in equation fl24j) . 



These scaling forms are fixed by the requirement that the density fluctuations obey the 
scaling law 

d 3 fc<|e o(M)l 2 > «r 4 (ii) 



with a similar law for the other components of 0^. We display our results for P p , X pU 
and P u in Figures |2|, |3| and (| along with fits to the following forms 

pP ^ = (i - (bz) + ( cz ) n y/ n (12) 

X pu (z) = t , r-^nr-r^i (13) 

v ; (1 - (ez) + (fz) m ) 2 / m v ; 

p u (z) = - (14) 

y } (1 - (hz) + (jz)P) om /P V } 

which are motivated by a requirement of white noise at large scales due to uncorre- 
cted initial conditions and good single parameter fits to P p ~ ^^-o.98±o.o6^ pu ^ 
(k£)- om±om and X pu ~ (^)" 2 - 0±ai at small scales. 

The errors indicate ensemble standard deviation. The other parameters are a = 1.18, 
b = 0.25, c = 0.24, n = 6, d = 0.23, e = 0.3, / = 0.5, m = 3, g = 0.07, h = 0.15, 
j = 0.18 and p = 1.59. We consistently get a peak for P p (z) and P u (z) at ~ 3 which 
corresponds to a physical wavelength of about 2£. However, it is difficult to be certain 
about the initial rise in the power spectrum because the error bars are large. In the 
notation of Magueijo et al the peak corresponds to an x c = krj of approximately 20. 
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Figure 3: Scaling fucnction for the equal time energy-momentum cross correlator as 
defined by equation (ffi). The solid line is the measured scaling function with 1 — a error 



bars from ensemble averaging. The dashed line is the fitted function equation (13). The 
dash-dot line is the predicted form from our model, equation (|36|). 
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Figure 4: Scaling fucnction for the momentum density as defined by equation (|T0| ) The 
solid line is the measured scaling function with 1 — a error bars from ensemble averaging. 
The dashed line is the fitted function equation flT4j). The dash-dot line is the predicted 
form from our model, the last line in equation (p6|). 
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As discussed above, the real importance of 0oo and 0oi in the context of perturbation 
theory is how they are correlated over time. We performed fc-space measurements of the 
two-time correlation functions 



C^(h 1 r ] ,r ] ') = V- 1 (Qoo(k,r ] )Q* 00 (k,r ] ')) 

C^(k;7 1 ,r ] ') = V- 1 (U(k,r ] )Q* 00 (k,r ] ')) 
C uu (k; V ,ri') = V- l (U(k, V )U*(k, rf)) 
During the scaling era they are well approximated by 



CP u 



C 



uu 



/£F 

1 

7T 



a 



T' 2 k(ri - rf)) 



P u (k£)P u {kZ)e 



■ ±r" 2 k 2 ( v - n ') 2 



l-T^e^-r]') 2 ) 



(15) 

(16) 
(17) 

(18) 

(19) 
(20) 



where we have used the scaling forms of the power spectra and cross-correlator and 
£ = £(v)j £' = £(v')- The imaginary component of these correlators is consistent with 
zero within the ensemble errors. The forms of these functions are motivated by a model 
which we describe in the next section. 

The values for T, T' and T" are given in Table |l|. A is approximately 35. These 
parameters were determined by minimisation of the %-squared for each realisation. The 
normal distribution of parameters obtained allows an estimate to be made of each pa- 
rameter and its 1 — o errors. We stress that we consider equations ([Of), fll9f) and (p0|) 
to be an approximation, although a good one for \r\ — vf\ < 8/c^ 1 . For comparison, the 



energy and momentum correlators fall to half their maximum value at \r) — rf 



3k~ 



Outside the range given there are two effects present in the measured correlators which 
invalidate the model. The first is small /c-dependent oscillations about zero as predicted 



by Turok ||22|| . The second is sharp peaks in the correlators at small scales which we take 
to be an effect of defining strings on the lattice. Figures |5|, ^| and [7| show the measured 



E c 


T 


a 






25 


0.21 ±0.05 


0.19 ±0.05 


0.42 ±0.05 


0.36 ±0.07 


45 


0.18 ±0.05 


0.16 ±0.05 


0.42 ±0.05 


0.42 ±0.06 



Table 1: Fitted parameters for the models in equations flilf) , ([19]) and ( p0|) 



functions at a time vf = 225 (which is within the scaling era of our simulations). 



Figure 5: Two-time correlation function (0 O o(&, 7/)0Q O (fc, 7/')) for rf = 225 




Figure 6: Two-time correlation function (U(k,r])QQ (k,r]')) for 7/ = 225 



4000 ^ 




Figure 7: Two-time correlation function (U(k,rj)U*(k,r)')) for 7/' = 225 
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From the gaussian fall-off over time we see that for a given mode with wavenumber 
k the network energy density decorrelates on a characteristic time scale which goes like 



l/TA;y(l — (kA)) and the momentum density varies as 1/T"k for E c = 25. 

4 Theoretical Model 

4.1 Power Spectra 



We can use a very simple model to account for the forms given in (|T2"D and (|H|), if not 
the precise details. From (H) we get 

6^(k) = J da (X»X U - x'^X'")^*^ (21) 

At any given time the power spectrum is 

(|0 OO (M)| 2 > = ( / da da' e 'Mx ((V? )-x( CT ',^ (22) 

We then make some assumptions about the statistics of the network. Firstly, that for 
each lag cr_ = a — a', the quantities Xi(a) — Xj(cr'), Xi(a) and X-(a) are gaussian random 
variables with zero mean. Secondly, that the distribution of strings is isotropic. Then 

(|6 00 (A;)| 2 > = fdada'e~ l e k2 ^- x ^ 

J (23) 



2 



da+da_ e -| fc2r (-) 



where we have introduced the function T(a — a') = ((X(cr) — X(<r')) 2 ) and changed 
variables to a_ = a — a' and a + = a + a'. We will also use t 2 = (X /2 ). The third 
assumption is that the network is described by a collection of randomly placed string 
segments of length £/i and total energy VC,^ 2 , where V is the volume of the simulation 
box. Then T = ?£ 2 (cx_/0 2 and L = \ J da + = V^ 2 . If we define z = aJb/£ then 

<|e o(M)| 2 > = Ljda.e-^ 2p ^~^ 2 

= H f\ dz e -m? z > (24) 

t J-\ 

V2VQ K N 
-erf(- 



& k£ K 2V6 ; 

Thus the scaling form of (@) emerges quite naturally. If we compare ( ^4[) with (||), we can 
write the predicted form of the scaling function P p (k^). The large and short wavelength 
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limits are 



P»(k£) = 



(k£)° if^«2 v / 6 
(k^)- 1 if^»2 v / 6 



(25) 



This function is plotted along with the measured and fitted model in Figure |2|. It should 
be noted that the normalisation is given by the model along with the measurement for t, 
and not introduced by hand. This model does better at predicting the measured power 
spectrum than a random walk model, in which the power on large scales goes like k~ 2 . 
Therefore, although each individual string is a random walk, they must be highly (anti) 
correlated on large scales. 

Similarly, the power spectrum for U becomes 



functions elsewhere ||19|1 , and we can use them to calculate the scaling function P u . 
The result is plotted in Figure |4j. Although for very large k£ this model predicts a 
(fc£) -1 behaviour, within the /c-range of our simulations the model agrees well with our 
measurements giving a fit to (k^)~ om , although the normalisation is not so impressive 
as for the energy- density. 

4.2 Two-time correlation functions 

In order to calculate the two-time correlation functions in this framework we need the 
two-time correlation functions T(a,o J ,rj,r} / ), V((T,a',T],T]') and n(cr, a', T], T]'). However, 
we have not yet measured these quantities and instead resort to crude modelling. If we 
consider small scales we may assume that each segment is moving at a velocity v 2 = (X 2 ) 
in a random direction orthogonal to the orientation of the segment so that 



(\U(k)\ 2 ) = (J dada'hkx'i^X^a^e^^-^'^) 
= Ljda. (|V(<t_) - §A; 2 nV_))e-^ 2r ^ 



(26) 



where V and II are the correlation functions 



V{a) = (X(a) ■ X(0)> 



(27) 




r(tr_, 77, r/) = Pal + v 2 (rj - r/) 2 



(29) 
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Using equations ( Pip and (|29| ) we get 

(900(^77)9^(^ 770) = (|e 00 (^r / )| 2 ) e -^ 2fc2 ^') 2 (30) 

However this model does not reflect the fact that between rj and rj' some energy is lost 
and the integrand is no longer independent of da + . We know that C pp (r], 77') = C pp {rj ', 77) 
and the most natural way of respecting this condition is to replace the power spectrum 
at 77 with the square root of the product of the power spectra at the two times. Hence 
the final result to be compared with equation (|18|) is 



(BooC*, 77)0^, 77')) = -±^y/pp(kt)Pr(k?) e- l e^-v'r {31) 

To model the other two-time correlators (U(k,T])U*(k,T]')) and (U(k,rj)QQ Q (k,r}')) we 
must be careful about the conservation of energy-momentum as loops are created and 
energy is lost from the long string network. If we assume that the loop production occurs 
evenly along the string and that we are in a scaling regime we can model in fc-space the 
rate of energy going into loops as 

A(M) = ^Me 00 (M) (32) 

The energy-momentum conservation equation becomes 

u(k,n) = ~e ao + ±e 00 {k,r,) (33) 



and hence, 



([/(M)9oo(M')> = / do " rf(j/ ^ + ^) e4fc2( ' 2(W)W) 

dada' (A _ \kv 2 ^))e-^ 2 ^ 2+p ^ (34) 

= (|e o(A;,r / )| 2 )(A_i^V) e -^ 2(W)2 

To ensure that we do not pick out a time 77 we express C pU (r),r)') 

(U(k, r))e* (k, rj 1 )) = (9oo(*, 77)9^, 77')) (—4= - \kv^J) (35) 



The time dependence of this function is plotted in Figure [8], and should be compared 
with Figure H. 
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Figure 8: Model of time dependence of (U(k, T])QQ (k, rf)) 




These plots assume a value of A measured from the simulations. We find it to be 
roughly constant at large at A ~ 0.32 (±0.03). Including loop production shifts where 
the cross-correlator goes to zero from 77 = 7/ to A — ^k 2 v 2 r]_^J (££')) = 0. The effect 
of energy loss through loop production on the momentum of the long string may also 
account for the form of the equal time cross-correlator. If we take 77 = 7/ in equation 
(|35l), which in terms of this model gives the correlation between the long string energy 
and the reaction momentum from loop production, we get 

X pU {K) = ^P P (kO (36) 

This neatly explains the {k^)~ 2 dependence in equation (|13|), although the measured 
normalisation is only 60% of that given by the independently measured A and equation 
(p6f). One reason for this discrepancy may be that our model assumes loop production 
occurs evenly along the string at a constant rate, whereas we observe loop production 
in more discrete bursts. The prediction from equation fl36|) is plotted in Figure ^, along 
with the measured scaling function. 

We find that the effect of the loop production term is insignificant for (U(k, rj)U*(k, rj')) 
and we will drop it in the following expression: 

(U(k,r])U*(k,T]')) = 1 fdada' |-A e -^V(W)W) 

k J or] dr]' 

= <|0oo(M)| 2 >y(l - jk 2 ( V - rWe-^W 

Again to ensure C uu (r),T]') = C uu (r]',r]), we use the product of the square root of the 
power spectra, so that 

(U(k,r])U*(k,r]')) = -^=^PP(kOP^')j(l - *#{fl-rff)e-t**^r (38) 

Equation ( pB]) incorrectly predicts U to have the same k£ dependence in its power spec- 
trum as the energy, with a relative normalisation given by v 2 /3. We expect a proper 
treatment in terms of the two-time on-string correlators for IT, V and T to remedy this 
problem. The time dependence of this function is plotted in Figure |Tl]. As one can verify 
by comparison with Figure [TT|, the simplest of assumptions has provided a reasonable 
approximation to the measured time dependence. 

In this section we have motivated the forms in equations (|18|) , ( |19D and (|20|) . From 
the model, and using the measurements t 1 = 0.637 ± 0.010 and v 2 = 0.363 ± 0.010, 
we can make predictions for the coherence parameters T, T' and T", and the relative 
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quantity 


model 


measured 


T 


0.35 ±0.02 


0.21 ±0.05 




0.35 ±0.02 


0.42 ±0.05 


'YV/ 


0.35 ±0.02 


0.36 ±0.07 


R 


0.35 ±0.02 


0.38 ±0.1 



Table 2: Parameters for the forms given in equations (0), ([19D and (|20|) as predicted by 
our model and measured from the simulations. Note a complication in comparing the 
values of T, due to a lattice effect in equation (]T8|) as discussed in the text. 



normalisation R = (Q 00 (k,r])Ql (k,ri')) /(U(k,rj)U*(k,r]')). The measured values are for 
E c = 25. The comparison for T is complicated by an apparent lattice effect accounted for 
in equation (|18|) by the length scale A. We assume that as A goes to zero and the form in 
equation (|T8"|) approaches that in equation (|31"D, allowing a comparison to be made. Not 
surprisingly, the model is not wholly satisfactory. We know that the momentum density 
does not have the same dependence as the energy density, which is predicted by (|38|). 
Also from we see that the strings decorrelate on a time scale ~ 1/k^J (1 — (kA)) and 
not k~ l . One might think that both features are understandable as it is well known that 
strings move faster on smaller scales. This gives faster decorrelation and greater power 
in the momentum density on small scales than the simple model would suggest. Having 
said this, we note that the scale A introduced in equation (|1~8"D is close to the lattice 
scale and is fairly constant throughout the simulation. Consequently, the departure from 
a k~ l coherence time may be due to the fact that the strings are defined at all times 
on the lattice. It is puzzling that this effect does not show up in the other correlators. 
Furthermore, the model does not account for the oscillations about zero in the two-time 
correlators. The oscillations are present in all three two-time correlators, but are most 
obvious in Figures ^] and [11]. As Turok has recently pointed out [22\, such oscillations 
should appear in the Fourier transforms of two-time spatial correlators because of a 
causality constraint which sets the correlator to zero for causally disconnected points. 



5 Conclusions 

We have demonstrated the scaling properties of the power spectra and cross correlator 
of two important energy-momentum quantities, the density 0oo and the velocity U, and 
given the large k behaviour. We have also studied the two-time correlators and measured 
the time coherence in the network. 
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We find that the energy and the momentum power spectra are peaked at around 
k£ ~ 3, where £ ~ 0.157?, thereafter decaying as an d (&£)~°' 66 respectively over 

the range of our simulations. The cross-correlator decays as {k^)~ 2 and is peaked at 
k£ ~ 2. 

For a mode of spatial frequency k, the two-time correlation functions of Goo an d U 
display correlations over time scales of rj c ~ 4.7k' 1 and rj c ~ 2.4/c _1 respectively. The 
time scale i] 2 is the variance in the Gaussian fall-off as a function of the time difference 
(see equations fljJD and ([20"D). We have presented simple models to explain the qualitative 
results, predicting the power spectra and cross correlator including normalisations, to a 
reasonable accuracy, and accounted for the features of the two-time correlation functions. 

The model describes the string network as a set of randomly placed segments of 
length where t — (1 — v 2 )^ , with random velocities. We assume that relevant 
quantities such as velocities and extensions between points with a given separation in a 
are gaussian random variables. We can then reduce ensemble averaging to the study of 
two-point correlations along the string, which we must model or measure. 

From the model we can show that the characteristic coherence time scale for a mode 
of spatial frequency k is r] c ~ 3/k. Our simulations give a similar result, although at very 
small scales r\ c decrease faster than k~ l . This we believe to be a lattice effect. 

There are potential pitfalls in taking the Minkowski space string network as a source 
for the fluid perturbation variables in a Friedmann model. For example, the energy 
conservation equation is modified to 

6>oo + -(©oo + Q)-kU = -A (39) 
a 

where a{rj) is the scale factor and 6 = Qu the trace of the spatial components of Q Mi ,. 
Since G is unconstrained by any conservation equation, its fluctuations could drive a 
fluctuating component in Gqo whose time scale would go like a/ a [2IJ. However, we find 



that (|©| 2 ) << (|Goo| 2 ) and believe the introduction of the a/a term to be a small effect. 

If we take the simulations at face value, the implications for the appearance of the 
Doppler peaks are not entirely clear cut. The coherence time is smaller than, but of 
the same order of magnitude as, the period of acoustic oscillations in the photon baryon 
fluid at decoupling, which is roughly 11 jk ||24j| . This is in turn smaller than the time 



at which the power in the energy and velocity sources peak, approximately 20 jk. The 
computations of Magueijo et al || of the Microwave Background angular power spectrum 
for various source models suggest small or absent secondary peaks. However, our string 
correlation functions can be used as realistic sources to settle the issue. 
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